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Abstract 

We investigate a family of distributions having a property of stability-under-addition, 
provided that the number v of added-up random variables in the random sum is also a 
\ random variable. We call the corresponding property a ^-stability and investigate the 

situation with the semigroup generated by the generating function of v is commutative. 

Using results from the theory of iterations of analytic functions, we show that 
the characteristic function of such a ^-stable distribution can be represented in terms 
of Chebyshev polynomials, and for the case of j/-normal distribution, the resulting 
characteristic function corresponds to the hyperbolic secant distribution. 

We discuss some specific properties of the class and present particular examples. 
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1 Introduction 



J> ■ In many applications of probability theory certain specific classes of distributions have 

become very useful, usually called "fat tailed" of "heavy tailed" distributions. The Stable 
distributions that originate from the Central Limit problem, are probably most popular 
£T) • among the heavy tailed distributions, however there is a wide collection of classes of 

00 . distributions, all related to Stable ones in many various ways, often these relations are not 

'. at all obvious. 

Besides, certain generalizations of stable distributions are known, using sums of random 
numbers of random variables (instead of sums with deterministic number of summands), 
see e.g. Gnedenko [3], Klebanov, Mania, Melamed [8], for the examples of such, including 
the so-called v-stable distributions, introduced independently by Klebanov and Rachev [9] 
■ and Bunge pp. 

In the present paper, we focuse on presenting further examples of strictly zv-stable 
random variables, that could be useful in practical applications, including applications in 
financial mathematics. 
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2 Definition of strictly //-stable r.v.'s, properties and examples 



In the present section, we give a general insight on strictly ^-stable distributions and 
describe some examples that have been mentioned in the literature before. 

2.1 Basic definitions 

Let X, X±, X2, ■ ■ ■ , X n , ... be a sequence of i.i.d. random variables, and let { u p , p £ A } 
be a family of some discrete r.v.'s taking values in the set of natural numbers N. Assume 
that this family does not depend on the sequence {Xj, j > 1}, and that, for A C (0, 1), 

E [u p ] = -, y p € A . (1) 
p 

Definition 2.1. We say that the r.v. X has a strictly v -stable distribution, if V p G A it 
holds that 

X^pV-j^X,, 
1=1 

where a € (0, 2] is called the index of stability. 

After this general definition, a narrower class is defined for a = 1/2. 

Definition 2.2. We call the r.v. X a strictly v-normal r.v., if EX = 0, EX 2 = 00, and 

the following holds: 

X±p^jrX 3 , Vp€A. 

Closely related to the stability property is the property of infinite divisibility, so we 
also give the following definition. 

Definition 2.3. X has a strictly v -infinitely divisible distribution, if for any p € A, there 
exists a r.v. Y^ p \ s.t. 

Up 

X = Y1 Y j P) ' with Y (P) ' ^ ' • • • ' y « P) ' ■ ■ ' bein9 iid r - v - ' s ■ 

A powerful tool for investigating distributions' properties is the generating function. We 
shall use the generating function of the r.v. v p denoting it by T , p (z) := E [z Vp ]. Moreover, 
we denote by A the semigroup generated by the family { V v , p £ A}, with the operation 
of the functions' composition. 

2.2 Summary of the known results 

With regards to the definitions above, the following results are known (see e.g. [7] for 
proofs and details). 

Theorem 2.1. For the family {V p , p £ A}, with E [up] = -, there exists a strictly 
v-normal distribution, iff the semigroup A is commutative. 
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Suppose that we have a commutative semigroup A. Then the following statements 
(that we refer to in the sequel as Properties) are known to be true (see [7j for proofs and 
details): 

1. The system 

<p(t)=V p (<p(pt)), VpeA. (2) 
of functional equations has a solution that satisfies the initial conditions 

<p(0) = l, ¥>'(0) = -l. (3) 

The solution is unique. In addition, there exists a distribution function (cdf) A(x) 
(with A(0) = 0) such that 

oo 

<p(t) = J e- tx dA{x). (4) 
o 

2. The characteristic function (ch.f.) of the strictly z/-normal distribution has the form 

f(t)=cp(at 2 ), a>0. (5) 



3. A ch.f. g(t) is a ch.f. of a ^-infinitely divisible r.v., iff there exists a chf h(t) of an 
infinitely divisible (in the usual sense) r.v., such that 

/(*) = ¥>(- In /»(*)). (6) 



The relation © allows obtaining explicit representations of ch.f. of strictly ^-stable 
distributions. Clearly, they are obtained through applying ([5]) to a ch.f. h(t) , provided 
that the r.v. corresponding to h(t) is strictly stable (in the usual sense). Moreover, note 
that the ch.f. ip(ait), a 6 K 1 , is the ch.f. of an analogue of the degenerate r.v., and that 
for the r.v. with such ch.f. the following analogue of the Law of Large Numbers exists. 

Theorem 2.2. Let X±, X2, ■ ■ ■ , X n , . . . be a sequence of iid random variables with the 
finite absolute value of the first moment, and { v p , p G A } a family of r.v. 's taking values 
in N, independent of the sequence {Xj,j = 1, 2, . . . }. Assume that E [u p ] = - and that the 
semigroup A is commutative. 

Up 

Then the series p Xj * s convergent is distribution, as p —> 0, and the limit of 

convergence is a r.v. having the ch.f. ip(ait). 

The proof of this theorem follows straightforwardly from the Property 1 outlined above 
and from the Transfer Theorem of Gnedenko, see e.g. [3]. 

In the following paragraph we discuss several particular examples of strictly z^-normal 
and strictly ^-stable distributions. 

2.3 Examples and the outline of the problem 

Example 1 . The usual stability. 
Assume the following setup : u p = - with probability 1, where p€A = |l,^,...,i,...}, 

and so V p {z) = z l / p . 
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Clearly, the corresponding semigroup A is commutative. 

oo 

Furthermore, <p(t) = exp{— t} = j e~ tx dA(x), where A(x) is a cdf with a single unit- 

o 

sized jump at x = 1. In this setup the strictly v- normal ch.f. is the ch.f. of the normal (in 
the usual sense) r.v. with the zero mean. 

Example 2 . The geometric summation scheme. 
Suppose, Up is the r.v. having a geometric distribution 

P{v P = k}=p(i- P ) k - 1 , k = l,2,... ,pe(0,l). 

Clearly, here E [u p ] = - , and V p (z) = 1 _^ z _ p -j z ,J)6 (0, 1) . It is quite straightforward 
to check that A is commutative. 

oo 

Moreover, a direct calculation gives tp(t) = = J e~ tx e~ x dx, i.e. A(x) is the cdf of 

o 

the exponential distribution. So that a ^-analogue of the strictly normal distribution is 
the Laplace distribution with the ch.f. f(t) = 1J ^ at i ■ 

Example 3 . Branching process scheme. 
Let V(z) be some generating function, with 'P'(l) = — > 1 (so that the introduced 
notation is po = 1/V'(1), with the condition po < 1). 

Consider now a family given by V 0n (z) = 7?°("-i) (-p(z)) , n = 1,2, . . . . Related to 

that is another family of the r.v.'s v v : V p {z) = V 0n (z) , p £ |^ , n = 1, 2, . . . j =: A . 

Clearly, the semigroup A coincides with the family {V p , p £ A} . The ch.f. ip(t) is a 
solution of the functional equation (p(t) = V (y>(pot)) . 

It can be noted that the content of the paper by Mallows and Shepp [11] is actually 
based on considering an example identical to the Example 3 above. Probably, neither the 
authors of that work nor its reviewers were familiar with the works by Klebanov and 
Rachev [9] and Bunge [I], which had dealt with exactly the same example a number of 
years earlier. 

Like mentioned in Introduction, in the present work we aim in widening the collection 
of examples that involve random summation with the commutative semigroup A. For that 
purpose, we address the description of pairs of certain commutative generating functions 
V and Q, i.e. the ones for which the balance equality V oQ = QoV holds, - but including 
only the case when there exists no such function % such that V = T-L ok and Q = T-L 0m for 
some k,m £ N (which would be exactly the case of the Example 3). 

In a general setting, the problem of describing all such commutative pairs of generating 
functions appears, unfortunately, far too involved to approach. However, certain special 
cases are rather straightforward for consideration. In order to approach the problem, we 
will use certain notions typical for the theory of iterations of analytic functions, that we 
outline in the separate section below. 

3 Theoretic justification via iterations of analytic functions 

Let V be a rational function with (deg) > 2. Denote by V 0n its nth iteration. The 
functions V and Q are called conjugates, if there exists a linear-fractional function R, such 
that V o R = R o Q . 

A subset E of the extended complex plane C is called completely invariant, if its 
complete inverse image V~ 1 {E) coincides with E. The maximal finite completely invariant 
set E(V) exists and is called the exceptional set of the function V . It is always the case 
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that card-E("P) < 2. Moreover, if card-E("P) = 1 then the function V is a conjugate to 
a polynomial, while for card.E('P) = 2 the function V is a conjugate to Q(z) = z n , n € 
Z\{0, 1} . Clearly, £(Q) = {0,oo}. 

If V is a rational function, then it is known (see e.g. (7j) that there is a finite number 
of open sets Fi, i = 1, . . . , r , which are Ze/t invariant by the operator "P and are such that 
(in the sequel, we will refer to the two points below as Conditions) 

r 

1. the union (J Fi is dense on the plane ; 

i=\ 

2. and V behaves regularly and in a unique way on each of the sets Fi . 

The latter means that the termini of the sequences of iterations generated by the points 
of Fj are either precisely the same set, which is then a finite cycle, or they are finite cycles 
of finite or annular shaped sets that are lying concentrically. In the first case the cycle is 
attracting, in the second one it is neutral. 

The sets Fj are the Fatou domains of V , and their union is the Fatou set F(P) of V . 

The complement of F(V) is the Julia set J(V) of V . Note that J{V) is either a 
nowhere dense set (that is, without interior points) and an uncountable set (of the same 
cardinality as the real numbers), or = C. Like F(V), J^P~) is left invariant by 

V, and on this set the iteration is rebelling, meaning that | V(z) — V(w) \ > \ z — w\ 
for all elements roina neighborhood of z (within J(V)). This means that V{z) behaves 
chaotically on the Julia set. Although there are points in the Julia set whose sequence of 
iterations is finite, there is only a countable number of such points (and they make up an 
infinitely small part of the Julia set). The sequences generated by points outside this set 
behave chaotically, a phenomenon called deterministic chaos. Let zq be a repelling fixed 
point of the function V, and let A = V'(zo). Define A : z — > Xz . Then there exists a 
unique solution of the Poincare equation 

FoA = VoF, F(0)=z Q , F'(0) = 1, 



that is meromorphic in C . 
Now let 

X{V) = F~ 

If for two functions V and Q we have V oQ = 

There are the two following possibilities: 
1. T(V) = C , in which case J(V) = C,. 



{J{V)). 

QoV, then they have the same function F. 



2. I(V) is nowhere dense and consists of analytic cuvrves. 

Fatou [2], and Julia [5] investigated the case. It turned out that is this case V and Q 
can be reduced by a conjugancy either to the form V{z) = z m and Q(z) = z n or to the 
form Viz) = T m (z) and Q(z) = T n (z) , where is the Chebyshev polynomial determined 
by the equation cos(kC) = Tfc(cosC)- 

4 Main results 

4.1 A new example 

Let us return to the study of i^-normal and z^-stable random variables. Recall that 
we deal with the family {u p , p G A} taking its values in N = {1,2,...}. As before, 
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we work with the generating function, V p (z) = E[z" p ],of v v . The important result 
that we stressed says the a strictly ^-normal (resp. strictly z^-stable) r.v. exist iff the 
semigroup A generated by {V p , p G A} is commutative. If V p , p G A , is a rational 
function (with deg < 2) satisfying Condition 2 of the above section, then either V p (z) 
is reduced to a form V p (z) = z x l p , p G {- , n = 1, 2, . . . }, and then we deal, in fact, 
with the classical (deterministic) summation scheme, or V p (z) is reduced to the form 
V p (z) = Ti/^p-(z) , p G { , n = 1, 2, . . . } . Clearly, the polynomial T m (z) is not a gener- 
ating function itself, however a function to which it is a conjugate, specifically the function 

W = - PG (I n = l,2,...| , (7) 



ri- 



'is indeed a generating function, - the fact that we prove below. Moreover, below we 
consider in some details a family of r.v.'s , p G , n = 1, 2, . . . } } that have generating 
functions of the form ([7]), and investigate the corresponding strictly ^-normal and strictly 
^-stable distributions. 

Lemma 4.1. Let P n (x) be a polynomial with degP n = k by the even powers of x, and 
whose zeros are all within the interval (—1, 1) . Let -P n .(l) — 1 an d polynomial's coefficient 
with with x n be positive. Then for any natural number k, the function 



V{x) 



2 k 



is a generating function. 
Proof. Represent P n {x) as 

n 

P n (x) = b + b\x H h b n x n = b n (x - aj) , 

j'=i 

where aj (j = 1, . . . ,n) are the zeros of the polynomial P n sorted in the order of ascen- 
dance. As P n is a polynomial by the even powers of x, then, if a.j is a zero of P n , then — Oj 
is also a zero of P n . Therefore, 

1 1 



1 



n/2 

1 " /2 1 

n/2 2 



^-Ilr^ (8) 



Obviously, 



x 



2 



a > 2 k=0 



E2 
a 3 



2fc 2fc+2 



X 
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is a series with positive (non- negative) coefficients, converging when |sc| < 1. From ([8j), 
it now follows that V(x) = x .\ is a series also convergent when \x\ < 1, having non- 

negative coefficients, and 'P(l) = 1 . Hence, V(x) is a generating function of some random 
variable. □ 

Corrolary 4.1. Let T n {x) be a Chebyshev polynomial of degree n. Then 

V(x) 



is a generating function of some r.v. taking values in N. 

Proof. When n is an even number, the result follows directly from Lemma and from 
the properties of Chebyshev polynomials. For odd n, consider the representation 
T n (x) = xP n -i(x), where P n -\{x) is a polynomial by the even degrees of x, satisfying the 
conditions of Lemma\4-1\ □ 



Let us now set A := {-^ , n = 1,2,... } . Consider the family of generating functions 

V P {z) = - — 1 — T1 , P eA. 



Clearly, V Pl o V p2 = V P2 o V P1 for all pi,p2 € A , due to the well known property 
of Chebyshev polynomials stating that T n (T m (x)) = T n . m (x). In other words, semigroup 
generated by the family {V p , p £ A} is commutative. It follows (see e.g. [7]) that there 
exists a solution to the system of equations 

v(t) = V p (tp(pt)) , peA, (9) 

satisfying initial conditions 

<p(0) = l, <f/(0) = -l, (10) 

and the solution is unique. 

Since T n {x) = cos(n • arccos x) = cosh(n • arccosh x) , the direct plugging gives that the 
function 

<p(t) = 1/cosh (\/2t) (11) 



satisfies the system ([9]), as well as the conditions (fTOjl . Hence, the function 



is actually a ch.f. of a strictly z/-normal r.v.. The ch.f. (|12p is, in fact, well known - it is 
the ch.f. of the hyperbolic secant distribution. Clearly, here a is the scale parameter. When 
a = 1, it is the case of the standard hyperbolic secant distribution, whose pdf has the form 

p( x ) = - S ech 



while the cdf is 

2 

F{x) = — arctan 

7T 



exp It) 



Furthermore, in order to obtain the expression for the ch.f. of strictly instable distributions, 
one just needs to apply the relation © to the strictly stable (in the usual sense) ch.f. h. 
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4.2 An interesting property 

Note that the function ip, as represented by (|lip . can be viewed somewhat interesting on 
its own, and so we shall address its properties and consider its cdf A{x) (which corresponds 
to <p(t) via (jD ) . 

Let Wi(t) and Wzit) , t > , be two independent Wiener processes. Consider a r.v. 

1 l 

f = J Wl{t)dt + j Wi(t)dt. (13) 
o o 

This r.v. is well studied, and it is known (see e.g. [12]) that its Laplace transform equals 
to 

E 



1 

cosh (V2t) ' 

which coincides with ip(t) as given by (jlip . 

Hence A{x) is the cdf of the r.v. £. On the other hand, as follows from Gnedenko's 
Transfer Theorem, 

A{x) = lim P { p v p < x } . 

p— >o 

Consequently, the following theorem is valid. 
Theorem 4.1. Let {v p < x , p £ A} be a family of r.v. 's having generating functions 

Vp{z) = —^— T - peA = \^,n = l,2,.. 

Then 

lim P {pv p < x} = P{^ < x} , 

p— >o 



where the r.v. £ is the one defined via {I 

Theorem 14 . 1 1 may be reformulated in the following way. 
Let 



1 

± n\ z ) fc=n 



k=0 

Then 

[n 2 x] 

lim N Pk{n) = P{£ < x}. 

n— >co — ' 

k=0 

On Figure 1, the plot of the X^fcLo' ^ s gi ven as a function of n starting with 
n = 2 until n = 50. We see that the functions attains the constant level rather quickly, 
and therefore it is possible to use the asymptotic result for n > 25. 

Corrolary 4.2. Let X be a r.v. having the standard hyperbolic secant distribution. Then 
its distribution can be represented in the form of a scale mixture of a normal distribution 
with zero mean and standard deviation -^/f , where £ is defined via U3\) . 
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10 20 30 40 



Figure 1: Plot of the X^L—o Pk{ n ) as the function of n = 2, . . . , 50 

°° 2 

To prove the above, one just needs to write the ch.f. of X in the form j e _t x dA(x) , and 

o 

note that e~ l x is actually the ch.f. of the standard Normal r.v. N(0, a 2 ) (a 2 = x), while 
A(x) is the cdf of £ . 

Note that there is a certain analogy between the representation A{x) as the cdf of the 
r.v. £ from (|13|) and the corresponding result in the scheme of the random summation with 
geometric distribution (see e.g. [7]). Specifically, considering the family {v p , pG (0,1)} 
having the geometric distribution P{f p = k} = p(l — , k = 1,2,... , the func- 

tion tp turns into 

oo 

j"^ = / e- tx dA 1 (x), 
o 

where A\(x) is the cdf of the exponential distribution, i.e. A\(x) = 1 — e~ x for x > and 
A\ = for x < . It can be checked that if 771 and 772 are two independent standard 
Normal r.v.'s, then A\ is a cdf of the r.v. £1 = r] 2 +77! > which is, in a way, related to (|13|) , 

4.3 Characterizations 

Let us now turn to the characterizations of the distribution of the r.v. (|13p and of the 
hyperbolic secant distribution. 

Theorem 4.2. Let X±, . . . , X n , . . . be a sequence of non-negative iid random variables, 
and Up, p G {-y , n = 2, . . . } , is a family of the r.v. 's having the generating function 
V p (z) = — — ; independent of the sequence {Xj , j > 1} . 

If, /or some /ixec? p € A , 

Xx^p^X,-, (14) 
(where "=" is the equality in distribution), then X\ has the distribution whose Laplace 
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transform is 



E ^ = ^M75)- a>0 - (15) 



Proof. The equality (|14p , in terms of the Laplace transform =Ee , can be repre- 
sented as 

= P p (¥(p*)) . (16) 

Clearly, the function 



cosh (\/ai) 

satisfies (jl6j) for any a > and, moreover, is analytic in the strip 1 1 \ < r ( r > ) . 

In the following, we use the results of the book by Kakosyan, Klebanov and Melamed 
[Bj. Example 1.3.2 of this book shows that { ^ a , a > } forms a strongly e-positive family, 
where e is a set of restrictions of Laplace transforms of probability distributions given in 
R+ on an interval [0, T] (0 < T < r) . 

Clearly, the operator A : f —> V p (f(pt)) on s is intensively monotone. 

The result follows from Theorem 1.1.1 of the above mentioned book (page 2). □ 

Theorem 4.3. Let X\, . . . ,X n , ... be a sequence of non-negative iid random variables, 
having a symmetric distribution, while {v p , p £ A} is the same family as in the previous 
Theorem. 

If, for some fixed p £ A , 

^i = p i/2 f;^, a?) 

then X\ has the hyperbolic secant distribution whose ch.f. is 

Proof. Quite analogous to the proof of the previous Theorem, with the difference that 
instead of Example 1.3.2, the use of the Example 1.3.1 from [6J is sufficient. □ 

5 Other examples 

There exist examples of the pairs of commutative functions, which are not rational. 
Here we refer to the two classes of such functions, the first of which was investigated by 
Melamed |lUj and the second appears at first in the present work. 

Example I . (See Melamed [10] for detailed study) 

Consider the family of generating functions 

pl/m z 

where p £ (0, 1) , and m is a fixed positive integer. Obviously, in the case m = 1, V p (z) 
reduces to the generating function of the geometric distribution, and has already been 
mentioned this case above. Hence, assume that m > 2. In that case, it is easy to check 
that 

" (4) = (TT^tF™ ' (20) 
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and therefore the ch.f. of the strictly i/-normal distribution (for the family {u p , p E A} 
having the generating function ()19|) ) has the form 



with a parameter a > . 



Example II . 

Consider the family of functions 

V v {z) = tt- , (21) 

where p£ { -y , n = 2, . . . } , and m > 1 (an integer) . 

Using a slightly modified version of the proof of Lemma 1, it is easy to check that 
V p (z) is a generating function of some r.v. v p , p £ , n = 2, . . . } for any fixed whole 
number m > 1 (surely, both V v and v v both depend on m, but we omit this dependence 
in the notation). 

The case m = 1 has already been considered above. For m > 2 analogous methods are 
applicable, and so will refer to the results only. Specifically, 

V>(t) = " 17^7 . ( 22 ) 

(cosh s/2m£) 

while the ch.f. of the corresponding strictly i/-normal distribution has the form 

f(t) = 1 — u ^ , (23) 

(cosh at) /m 

where a > . 

Note that in the case m = 2, we have the following expressions for the distributions 
whose Laplace transforms are (f20|) and ([22]) . 
For m = 2, the formula (|20p gives 



% /TT2t 

This function is the Laplace transform of the distribution of the r.v. X 2 , with X being 
the standard Normal r.v. 

In a similar way, (j22[) gives for m = 2 



\/ cosh V4i 



l 

This function is the Laplace transform of the distribution of the r.v. I = J X 2 (t)dt , where 

o 

X(t) is the standard Wiener process. 
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